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Modeling the total and polarized emission in evolving galaxies: "spotty" 
magnetic structures 
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Future radio observations with the Square Kilometre Array (SKA) and its precursors will be sensitive to trace spiral 
galaxies and their magnetic field configurations up to redshift z ^ 3. We suggest an evolutionary model for the magnetic 
configuration in star-forming disk galaxies and simulate the magnetic field distribution, the total and polarized synchrotron 
emission, and the Faraday rotation measures for disk galaxies at 2; ^ 3. Since details of dynamo action in young galaxies 
are quite uncertain, we model the dynamo action heuristically relying only on well-established ideas of the form and 
evolution of magnetic fields produced by the mean-field dynamo in a thin disk. We assume a small-scale seed field which 
is then amplified by the small-scale turbulent dynamo up to energy equipartition with kinetic energy of turbulence. The 
large-scale galactic dynamo starts from seed fields of 100 pc and an averaged regular field strength of 0.02 /xG, which then 
evolves to a "spotty" magnetic field configuration in about 0.8 Gyr with scales of about one kpc and an averaged regular 
field strength of 0.6 /iG. The evolution of these magnetic spots is simulated under the influence of star formation, dynamo 
action, stretching by differential rotation of the disk, and turbulent diffusion. The evolution of the regular magnetic field 
in a disk of a spiral galaxy, as well as the expected total intensity, linear polarization and Faraday rotation are simulated 
in the rest frame of a galaxy at 5 GHz and 150 MHz and in the rest frame of the observer at 150 MHz. We present 
the corresponding maps for several epochs after disk formation. Dynamo theory predicts the generation of large-scale 
coherent field patterns ("modes"). The timescale of this process is comparable to that of the galaxy age. Many galaxies 
are expected not to host fully coherent fields at the present epoch, especially those which suffered from major mergers or 
interactions with other galaxies. A comparison of our predictions with existing observations of spiral galaxies is given and 
discussed. 
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^ 1 Introduction 



Magnetic fields of nearby galaxies demonstrate several mag- 
netic configurations and related configurations of polarized 
radio emission including a "ring" of polarized emission in 
M31, M 33 Ir^batabaei et al.ll2QQ8l) and magnetic arms in 
NGC 6946 (lBeckll2QQ5h . a" particular magnetic configura- 
tion in a galaxy is believed to be determined by specific fea- 
tures of the excitation mechanism. In other words, magnetic 
configurations trace important details of galactic MHD. For 
example, the magnetic configurations of the barred galax- 
ies NGC 1097 and NGC 1365 trace gas flows presumabl y 
feeding the central active galactic nuclei (B eck et al ] l2QQ5h . 

Future radio observations with the Square Kilometre Ar- 
ray (SKA) are expected to give rich information concern- 
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ing magnetic configurations in the disk of the first g alaxies 
with redshifts up to at least z = 3 dArshakian et al.l l2009). 
The aim of this paper is to suggest a plausible evolution- 
ary model for the magnetic configuration in early galaxies 
and to simulate the magnetic field distribution, and polar- 
ized synchrotron emission, and the Faraday rotation mea- 
sures for disk galaxies aX z ^ 3 (the polarization survey in 
the frame of the SKA Design Studies) which can be tested 
by future observations with the SKA and its precursors. 



On first sight, the answer on the question under discus- 
sion looks quite straightforward. The mean-field dynamo is 
an amplification and ordering mechanism which excites one 
or few modes of regular magnetic fields. The mode with 
maximal growth rate should dominate very soon, indepen- 
dent of the configuration of the seed magnetic field, and 
determine together with the distribution of relativistic and 
thermal electrons the pattern of polarized emission. Later, 
the magnetic field becomes dynamically important and the 
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dynamo becomes nonlinear, while the main shape of the 
configuration survives. 

A closer look, however, shows that it is much more com- 
plex. The first point is that galactic disks are thin and the 
time required for a dynamo to establish the magnetic field 
pattern perpendicular to the galactic plane, tw, is much shorter 
than the corresponding time required to el aborate the con- 
figura tion in radial extent of the disk, tr (lArshakian et al. 



2QQ9h . (We use a cylindric coordinate system r, 0, w con- 
nected with the equatorial plane of the disk.) In particular. 



the diffusion time scales with distance as 



If we as- 



sume that tr and are proportional to the corresponding 
diffusion times, then tr/t^ = {R/h^ ~ 2 x 10^ for a disk 
aspect ratio h/R^ 2 x 10~^ (here h and R are the vertical 
and radial scalelengths of the ionized gas). 

The second point is that if the seed field was not ordered 
at the overall galactic scale, the growing regular magnetic 
field can preserve alternating directions in l arge r egions for 
a long time. As suggested by IPoezd et al.l ( 1993 ) (see also 
Beck et aP 1 19941) . a suitable seed magnetic field in galac- 
tic disks can be produced by the fluctuation dynamo. It is 
then natural to expect that a large-sca le magnetic s tructure 
with rev ersals along radius (Kuzmaikin et al.l 1 1985k and/or 
azimuth (iBykov et al.lll997 ) can persist for a period compa- 
rable to the galactic lifetime. 

The third point is that tuj has to be shorter than the galac- 
tic lifetime. However, tr may be comparable to the galactic 
age and the magnetic field becomes dynamically important 
before the dynamo develops a unique magnetic configura- 
tion over the whole galaxy. Then the spotty magnetic field 
structure evolves nonlinearly and can remain complicated 
for a rather long time. Observations of magnetic fields in 
nearby galaxies confirm that their magnetic configurations 
are more c omplicated t han just a simple leading dynamo 
eigenmode (lBeckll2QQ5h . 

In this paper, we try to address these points in a simpli- 
fied way. Summarizing, observations of magnetic configu- 
rations and their evolution with redshift z can be a tracer 
for the structure of the seed magnetic field and its origin. 
Obviously, events like encounters, mergers, etc. make this 
picture much more rich and also more complicated but this 
is beyond the scope of this paper. 

In Sect. [21 we present the physical bases of the dynamo- 
evolution model. The model of the evolution of magnetic 
fields in disk galaxies is described in Sect. [3] The link be- 
tween star-formation rate and total magnetic field of a galaxy 
is presented in Sect.lH and modeling of the regular magnetic 
fields, total intensity, polarization and Faraday rotations of 
evolving galaxies are presented in Sect. |5] Discussion and 
conclusions are given in Sect. [6l 

2 Physical picture of evolving regular 
magnetic fields 



A three-phase model of the evolut ion of regular niagneti c 
fields in galaxies was developed by lArshakian et al.l (|2009|) . 



The dynamo theory, which successfully describes the polar- 
ized synchrotron emission and Faraday rotation in nearby 
galaxies, was used to derive the timescales of amplification 
and ordering of magnetic fields on small and large scales 
in disk and quasi- spherical galaxies. This allows predicting 
the generation and evolution of magnetic fields in young 
galaxies at high redshifts. We refer to disk galaxies that 
rotate differentially and the turbulence is driven by super- 
novae, and we adopted the standard values of turbulence 
velocity v = 10 km s~^ and a basic scale of turbulence of 
/ = 100 parsecs. Results from simulations of hierarchical 
structure formation cosmology provide a tool to develop an 
evolutionary model of regular magnetic fields coupled with 
galaxy formation and evolution. 

In the first phase (z ^ 20), the seed magnetic fields 
of order :^ 10~^^ G (lArshakian et al.lEQQ9[ and references 
therein) were generated in dark matter galactic halos by 
the Biermann battery mechanism or the Weibel instability 
(Weibel 1957; Schlickeiser & Shukla 2003; Lazar et al. 2009). 
Weak seed magnetic fields can be exponentially amplified 
by th e small-scale dynamo during the fonn ation of the first 
stars dSchleicher et al.l2Q10l:ISur et al.l2Qlol) and by the Bier- 
mann mechanisni in su pernova explosions of the first stars 
( Hanayama et al ] l2005l) . 

The key feature of this model is that, in the second phase, 
the accreted gas from the intergalactic medium and merg- 
ers of small-mass dark-matter halos generated turbulence 
on scales comparable with the size of the infalling gas in 
the halo of a protogalaxy. The small-scale dynamo driven 
by turbulent, shock-heated gas in the halo of a protogalaxy 
can rapidly amplify the seed field from ~ 10~^^ G to the 
equilibrium level, ~ 10~^ G, in a relatively short time, ~ 
4 X 10^ years. This pha se was for the first time proposed by 
lArshakian et al.l (|2009|) and it may happen at z ~ 10 — 20 
or even earlier. 

In the third phase, the extended thin disk was formed at 
redshifts from z < 10. The formation of gas-rich disks can 
be triggered by the dissipation of the protogalactic halo or 
by major merger events. A weak, large-scale magnetic-field 
component was generated in the disk from small-scale mag- 
netic fields of the halo (amplified to the equilibrium level in 
the second phase) and further amplified to the equipartition 
level by the mean-field dynamo within few billion years. It 
was assumed that the small-scale galactic dynamo is acting 
together with the large-scale dynamo and that it saturates at 
the equiparti tion with the kinetic energy of interstellar tur- 
bulence (e.g. ISubramanianll 1998k . The evolution of the co- 
herence scale of the regular large-scale magnetic fields took 
a longer time: full coherence was reached in dwarf galaxies 
within ~ 6 Gyr and in disk galaxies within ~ 8 Gyr. Giant 
galaxies need more than 15 Gyr to develop a fully coherent 
field. 

The above picture is not the only possibility. The seed 
magnetic fields could as well be produced by the turbulent 
dynamo in an already formed galaxy as soon as turbulence 
driven by star formation develops (more precisely, in 10^- 
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10^ years after the first star formation burst). Then the first 
two stages are unimportant, but the form of the seed mag- 
netic field for the large-scale dynamo action remains the 
same as in the above picture (but delayed). Therefore, in- 
dependent of the poorly known details of galactic formation 
and their early evolution, we can more or less confidently 
suggest a plausible and robust distribution of the seed mag- 
netic field (see Sect. 13. II) . 

The dynamo operates in the disk and is able to amplify 
the amplitude of t he regular field and to o rder the field on 
timescales given in lArshakian et al.l ( 2009 ). 



3 Model of evolving magnetic fields in 
galaxies 

3.1 Origin of seed magnetic fields in first galaxies: 
magnetic spots 

The large-scale galactic dynamo starts from a configuration 
disordered on the scales of the disk radius and inevitably 
gives a "spotty" distribution of magnetic fields over the galac- 
tic disk. Magnetic fields in nearby spots have opposite di- 
rections (clockwise and counterclockwise). It looks plausi- 
ble that the spotty nature of the magnetic field (see Fig. 3) 
can be isolated observationally, provided the resolution and 
sensitivity allow observations in galaxies with suitable red- 
shifts. 

There is however an alternative viewpoint. The early 
Universe could contain magnetic fields. This field could be 
homogeneous and strong enough to serve as a seed field for 
the galactic dynamo, while weak enough to keep the Uni- 
verse practically homogeneous and isotropic. Then the seed 
magnetic field is coherent on the scales much larger than 
the galactic size, and one can hope in principle to avoid a 
phase with "spotty" magnetic structure and obtain a mag- 
netic field ordered on the scale of the whole galaxy at all 
stages of galactic evolution. This scenario looks less attrac- 
tive, but it cannot be rejected by the available observations 
and theoretical ideas. Moreover, the existence of a homoge- 
neous cosmological magnetic field gives an important mes- 
sage for cosmology. Causal mechanisms for magnetic field 
generation in the early Universe, based on standard practical 
physics, result in magnetic fields with contemporary scales 
which are negligible compared to galactic scales. This means 
that the desired primordial homogeneous magnetic field wo- 
uld have an a-causal nature, i.e. the scale was larger than the 
horizon at the time of its creation. 

It looks potentially possible to distinguish between the 
above scenarios by future observations with the SKA and 
hence support or reject the option of an a-causal primor- 
dial magnetic field. In case of the observationally confir- 
mation of the a-causal idea, primordial magnetic field will 
be supported observationally, the concep t of standard parti- 
cle physics have to be reco nsidered (e.g. iGiovanninilEooi ; 
Semikoz & SokolofflEoQSh . 



Obviously, both scenarios require to be elaborated. To 
be specific, we present here the results concerning the first, 
most probable scenario while the evolution of the homoge- 
neous primordial seed field inside a galaxy remains for the 
time being out of the scope of this paper. 



3.2 Initial magnetic spots: configuration and evolution 
of the radial and azimuthal ordering scales 

We start the description of the model from phase 3 when the 
disk is formed (t = 0). Assuming that the magnetic field 
in the disk is produced by the turbulent dynamo, we start 
with an rms turbulent magnetic field 6 = 6.2 /iG assuming 
equipartition with kinetic ener gy of turbulence, and whose 
scale radius is Iq = 100 pc ( Arshakian et al ] l2009h . Then 
the average amplitude of the regular field over the scale Iq 
at t = is given by 



Bo - bN- 



-1/2 



0.02/iG, 



(1) 



where A^i = ?>hR^/{2ll) 75, 000 is the number of turbu- 
lent cells in the volume of the disk for a galaxy with Gaus- 
sian radius scale of R = lOkpc and Gaussian height scale 
h ^ 500 pc. We assume that the regular field in initial mag- 
netic spots is oriented randomly in the plane of the disk - the 
vectors Bq change randomly from one spot of scale radius 
Iq to another. 

The regular magnetic field is spatially ordered on a co- 
herence scale of Iq = 100 pc. The radial size of a magnetic 
spot grows with time due to dynamo action. An approximate 
estimate of the ordering scale of the spot in radial direction 
is given by 



lr{t) = lo^Vrt, 



(2) 



where t is the age of a galaxy from the disk formation epoch, 
and Vr is the speed of the magnetic field pro pagation of the 
spot in radial direction. According to M oss et al.l (Il998k . it 
depends on the radial turbulent diffusivity ((3r = lov/3) and 
the dynamo growth rate (F = ftlo / h): 



Vr 



(3) 



where ft and h are the angular rotation and scale height of 
the disk, and ^ is th e turbulence velocity of the gas (see 
lArshakian et al.ll2QQ9[ and references therein). 

The turbulent cells as well as the magnetic spot s are 



stretched in azimuthal direction by differential rotation (iRuzmaikin et al 
Il988k . i.e. r and r + Ir move with different rotation rates 
ft{r) and ^}{r -\- Ir). However, the effect of turbulent diffu- 
sion in course of dynamo action works against stretching of 
a spot in azimuthal direction in the way that the pitch angle 
(p) of magnetic lines remains constant. We suppose that the 
turbulent diffusivity prevents further stretching of the spot if 
the ratio 1(f) /Ir becomes larger then tan~^ p, i.e. if the spot 



www.an-journal.org 



© 2006 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



4 



T.G. Arshakian et al.: Modeling the magnetic fields 



diagonal becomes more inclined than the magnetic line. As- 
suming that ^}{r) — rt{r -\- Ir) ~ Vt{r)lr/r we arrive to the 
following ordering scale in azimuthal direction 



at the age t, the amplitudes of the turbulent and regular fields 
can be measured from the set of equations 

5totW' = 62+5(t)^ (6) 



lr{t) 1 + 



h/B(t) = 0Viexp(-t/r), 



(7) 



lr{t) tan ^ p, if ( 1 



> tan p. 



We stress that the relations for radial and azimuthal or- 
dering scales are derived from the assumptions of our qual- 
itative model and that the ordering scales are likely to be 
model dependent. A specification and improvement of scal- 
ing laws on the basis of numerical simulation of a particular 
model of the galactic dynamo is highly desirable. 

The timescales of ordering the field on scales of 1 kpc 
in vertical and radial directions are comparable 0.8 Gyr; 
Eqs. (10) and (11) in Arshakian et al. 2009). We start the 
modeling of magnetic spots at the age of 0.8 Gyr from the 
scale radius of L = 0.5 kpc (scale length of 1 kpc) which 
has a size of the full height of the galaxy disk, 1 kpc. At 
this epoch, only a fraction of spots will contribute to the 
dynamo-driven large-scale magnetic field. Thus the num- 
ber of useful spots of a scale radius L = 0.5 kpc are those 
which have a magnetic pitch angle of the growing large- 
scale magnetic field, po ± Sp, where Sp is some tolerance, 
and the directions of small-scale magnetic fields in cells b 
are uniformly distributed in pitch angle p. The fraction of 
useful spots is 2Sp/27T, so that the number of useful spots is 



(4^here t* = h/{Ql) is t he dynamo timescale for amplifica- 
tion of the regular field (lArshakian et al ]|2009[ Eq.(9)),and 
y/N^ = b/B{t = 0) (see Eq. ([T])) is the number of turbulent 
cells in the disk of a galaxy at t = and depends on the ra- 
tio of geometrical volumes of the disk and the turbulent cell. 
At the epoch of a disk formation, the ratio b/B{t = 0) ^ 
300 for a disk galaxy with R = 10 kpc, h = 500 pc and 
/ = 100 pc, and b/B{t ^ 13.5 Gyr) = 1.7 for a contempo- 
rary Milky- Way type galaxy (MW) at the present epoch. If 
the age of a galaxy is smaller than the time when the regu- 
lar field reache s the equilibrium level, tgqi = t*/n(2/v^]Vi) 
(see Eq. (17) in lArshakian et al.ir2009|) . then we derive 



Btot{t) 



7m 



and 



B{t) = Btot{t)- 



^t/t* 

7m 



1 



ot/t* 



lft> teqh we assume that b/B{t) 
2 



^ 2 and hence 



V5 



TT 



(5) 



Btot{t), 

Btotjt) 
v/5 • 



(8) 



(9) 



(10) 



(11) 



3.4 Topology of the magnetic field 

The spatial structure of the total galactic magnetic field at 
any age is modeled as a superposition of a regular field B(t) 
in the disk and a random field b which describes the contri- 
bution of the large-scale galactic turbulence. 

To describe the regular magnetic field B we apply the 
following parametrization. 



Br (r, ^, w) = B{r^ ^, t) exp 



B^{r,(j),w) = 5(r, 0,t)exp 



where the factor 2 allows for the fact that both B and —B 
can seed the large-scale dynamo if B has a pitch angle of 
p ± 6p, and N2 — (R/L)'^ = 400 is the number of all spots 
in the disk. Experience with fitting azimuthal modes to the 
observed magnetic fields indicates that the magnetic pitch 
angle has, typically, a tolerance range Sp 5°. This esti- 
mate can be used, rather arbitrarily, in our case. Taking these 
numbers, we get n 20. Thus, about 20 magnetic spots 
have the pitch angles within 15° and 25° and constitute a 
regular field of B{t) averaged over a radial scale of 1 kpc . 
According to the evolutionary model Arshakian et al.l2009 ) 
the regular field of B{t) = 0.6 jiG will be achieved in t ^ 
0.8 Gyr (or at 2: 4.7) if the thin disk formed at z = 10 
(see Sect. Hand Table©. 

3.3 Evolutionary equations for turbulent and regular 
fields 

During the evolution of an undisturbed (by major merg- 
ers) disk galaxy with a slowly changing SFR, the turbulent 
field (6) remains almost unchanged, while the regular field 
is amplified to the equilibrium level 5(teqi) (where teqi is 
the time when the strength of the regular field reaches the 

equilibrium level) in a few Gyr and remains at this level field at the age t (Eqs. (1911 II) ), lr{t) and /^(t) are the ra- 
(lArshakian et al .l2009h . Given the total magnetic field strength dial and azimuthal ordering scales of the regular field at t 



sinp, (12) 



cosp.(13) 



The amplitude of the magnetic field in the middle plane 

21 



B(r,(t),t) = 5(t)exp - 



{(j) - (j)Q)rQ 



exp 



ro 



lr(t) 











exp 









tanh 



2' 



(14) 



where ro , 0o are coordinates of the position of the magnetic- 
field spot in the disk, B{t) is the amplitude of the regular 
field at the age t (Eqs. 
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(Eqs. © and Q), rmf, hmf are the Gaussian radius and 
vertical scales of the magnetic galactic disk, and p is the 
pitch angle of the magnetic field lines. We exploit the Carte- 
sian coordinates (x, w) or cylindrical coordinates (r, ^, w; 
X = rcos^, y = rsin^). The coordinates of the galaxy 
projected on the sky are (x^y' = ycosi), where i is the 
inclination angle of the galaxy. 

The first two exponential terms in Eq. (O describe the 
radial and azimuthal profiles of the spot, respectively, the 
forth term means that there is a cut-off in spot distribution 
nearby the outer border of the galaxy, and the last one stands 
to make the field amplitude vanishing just in the galactic 
center to make it smooth at this exceptional point. 

We found the vertical magnetic field component Bu, in 
the disk from the solenoidility condition 

dB,, 1 rOrBr dBs 



dw 



r \ dr 



(15) 



The turbulent magnetic field is considered as a diver- 
gence free, random fluctuating field with the given energy 
spectra ^ and shape function 11: 



(16) 



6(r, (p^w) = b ^(r, w) n(r, (/), w) 



where b is the strength of the turbulent field (Eqs. (ISIlOt ) 
and 



n(r, (/), w) = exp 



fe)1 hte) 



(17) 



where the spectral property of the random function ^ are 
specified as follows: 

{k/k^Y , k > ko 
{k/kof, k<ko. 



(18) 



For the sake of definiteness, we ad opt a = —5/3 Kol - 
mogorov scaling), /3 = 2, /cq = ( Stepanov et al.l2QQ8 ). 



3.5 Other assumptions of the model 

The assumptions of t he model for the magne tic field evo- 
lution is described in ^ Arshakian et al.l ( 2009 ). An average 



gas density of a protogalactic halo of « 6 x 10~^^ gcm~^ 
within its virial radius of about 600 pc was assumed at z = 
15. In the case of equipartition between magnetic and ki- 
netic energy of turbulence of lOkms"^, the gas densities 
of nearby disk galaxies varies between 2 x 10~^^ gcm~^ 
and 18 x 10~^^ gcm~^, corresponding to the range of ~ 
(5 — 15) /iG for the strength of the random magnetic field 
(IBeck et al.ll2005|) . For simplicity, we assume a gas density 



of p = 10~^^ gcm~^ that is unchanged during the evolu- 
tion. 

For a singular isothermal sphere, the rotation curve of 
the halo is given by Vc — 10i^(2:)r2oo, where is the cir- 
cular velocity of the dark halo at a limiting radius of r2oo 
within which the mean mass density is 200 times the back- 
ground density at the redshift z , and H{z) is the Hubble 
constant at redshift z (IMo et al.| [l998). If the gravitational 
effects of the disk itself are neglected (the thin-disk rotation 



(19) 



follows the rotation of the halo) then the angular rotation 
curve of the disk is 

K ^ lQij^(^)r2oo 
r r 

In the disk-halo model ("Mo et al.l ll99"8l) r2oo = Tc ^2/0.04 
and Vc = 350H{z)rc where Tc is the scale radius of the 
total cold-gas component (molecular and neutral hydrogen; 
H2+HI) which has an exponential distribution in the disk. 
As a first approximation, we assume that H{z) = H{0) = 
70 km s~^ Mpc~^ is constant, and that the angular rota- 
tion of the Milky-Way type galaxy has a unique curve at 
all epochs, 

24 5r 

l][km sec-^kpc"^] = — ^. (20) 
r 

According to observations, the exponential scale radius 
of the thermal emission (proportional to n^) is about the 
same as the exponential scale radius of the synchrotron in- 
tensity (e.g. I Walsh et al.lE002|) . so that the exponential scale 
radius of the thermal electron density Ue is related to that 
of the synchrotron intensity /gyn as In^ = 2/syn- We assume 
that the exponential scale radius of the total cold gas in a 
disk galaxy is /c(H2+HI) = 5kpc, which is similar to /gyn 
and the scale radius of the thermal gas /th, = ^syn = 
/th = 5 kpc. The scale radius of the electron density ng is 

= 10 kpc (/n2 = kh)- In this paper, we use the Gaussian 
scale radius which is r = /Vln 2 for distributions with the 
same half-power width. The relations between the Gaussian 
scale radii of the synchrotron intensity rgyn, cosmic ray rcr, 
and the magnetic field rmf (which is assumed to be the scale 
radius of all field components: total, regular and turbulent) 
follow from energy equipartition between magnetic fields 
and cosmic rays (^tot ^ ^cr) 

rc r^yn = r^i 4.2 kpc, 



rmf = 2 r. 



syn 



syn 



5.9 kpc, 
- 5.9 kpc, 
B.Skpc. 



(21) 



The exponential scale height of the Reynol ds layer of 

ionized gas (Ue) is about 1 kpc in the Milky Way ( Savage & Wakkei 
2009h . As this is hard to measure in external galaxies, we 



assume this height to be universal. This value is about one 
tenth of the exponential scale radius, and we may assume 
the same ratio for the Gaussian scales dStepanov et al.l2008|) . 
The exponential scale height of synchrotron emission is abou t 
1.8 kpc for the thick disk and 300 pc for the thin disk (IKrause 
2009). Dynamo action occurs in both components. The av- 
erage value of about 0.5 kpc is again about one tenth of the 
exponential scale radius of synchrotron emission. Hence, 
Gaussian scale heights can be assumed to be equal to one 
tenth of the scale radius for all components. 

We further assume for evolving MW-type disk galaxies 
(^ 10^^ Mq) that (a) the turbulence scale and turbulence 
velocity were driven by SN explosions and had values close 
to those observed in present-day disk galaxies such as the 
MW, / = 100 pc and v = lOkms"^; (b) the scale height of 
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the disk at present epoch is = 0.5 kpc and the scale length 
is R = 10 kpc, the disk is thin (h/R = 0.05) and the ratio 
h/R is constant at any epoch; (c) the scale length is constant 
over the cosmological epoch (no-downsizing); (d) the pitch 
angle of the regular field is unchanged during the evolution 
andp = 20°. 

4 Link between star-formation rate and 
amplitude of the total magnetic field 

In our paper, we do not introduce a star-formation rate (SFR) 
as an additional parameter of the model but rather link it 
to the amplitude of the total magnetic field. This can be 
used to scale the total and polarized intensities modeled in 
Sect. [5] The SFR is an important parameter which can be de- 
termined from ultraviolet (UV), Ha emission line, infrared 
(IR) and radio surveys. The UV and optical SFR indicators 
are strongly affected by dust, and therefore the IR and radio 
luminosi ties ha ve been used to estimate an unbiased SFR 
(see|M3 I2003L and references therein). 



Star formation rate 



IR - Radio 



Radio luminosity at 1.4 GHi 



of the tight radio-IR correlation, this implies that the radio 
synchrotron emission is also suppressed in low-luminosity 
galaxies, compared to brighter galaxies which already has 
been ob served in nearby spiral galaxies (Chvzv et al. 20071 



Dumke et al. 2000) and in dwarf galaxies ( Klein 1991). This 



is taken into account in deriving the relation between the av- 
erage SFR (in MqVT" ^)and the rest- frame radio luminos- 
ity at 1.4 GHz rsee lBd][200l Eq. (6)): 



SFR(radio) : 



5.52 X lO-^^Li.4, if Li.4 > Lc 



5.52 X lO-^^Li.4 



(22) 



if Li.4 < Lc 



^ 0.1 + 0.9(Li.4/Lc)0-3 
where = 6.4 x 10^^ W Hz-^. 

From the observed radio intensity of a galaxy one can 
estimate the radio luminosity and hence the star formation 
rate from Eq. (|22]) . On the other hand, the radio intensity 
can be used to estimate the total magnetic field strength 
( Beck & Kr ause 2005) under the assumption of the equipar- 
tition between the energies of the CRs and the total mag- 
netic field. In turn, the regular and turbulent magnetic-field 
amplitudes can be estimated from the total field strength 
(Eqs. ([SlfTTt) using the timescale of amplific ation of the reg- 



ular fi elds of the dynamo-evolution model ( Arshakian et al. 



l2009h . The scheme presented in Fig.[T]demonstrates the link 
between the SFR and the amplitudes of the total magnetic 
field, regular and turbulent fields. 



Redshift 



Radio intensity at 1.4 GHz 



Equipartition energy 



Total nnagnetic field (S^J 



Dynamo evolution 



Regular and turbulent fields 
(S and b) 



Fig. 1 The diagram shows the link between SFR and ma^ 
netic field strength in a star-forming galaxy. 



Here, we derive the link betwe en the SFR and the am- 
plitude of the total magnetic field, lid] (l2003h used radio 
and IR luminosities of star-forming galaxies to calibrate the 
radio-derived SFRs. He showed that IR emission is a good 
indicator of the star formation in luminous galaxies, while 
in galaxies fainter by two orders of magnitude the IR emis- 
sion traces a small fraction of the star formation. Because 




Fig. 2 Comparison of SFRs of 14 nearby galaxies esti- 
mated from the radio emission at 1.4 GHz and the IR or 
Ha emission. SFRs of nine galaxies (circles) are taken from 
Young et al.l (Il989h and five galaxies (triangles) from other 
papers (see Table [TJ. The dashed line represent the equality 
and the full line is the best linear fit to the data. 



To test this scheme, we selected a sample of 14 nearby 
spiral galaxies with different inclinations and levels of star 
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Table 1 Radio luminosities, total field strengths and star-formation rates (SFR) in nearby spiral galaxies 



Source 


Type 


i 


D 


S4.85 




i^tot 


SFR(radio) 


SFR(IR, Ua) 


Ref. 






n 


(Mpc) 


(mJy) 


(W Hz-^) 


/iG 


Moyr-i 


Moyr-i 




IC 342 


SAB(rs)cd 


25 


3.1 


856 


2.5 X 10^' 


10 


1.6 


1.2 


1 


M31 


SA(s)b 


76 


0.7 


1863 


2.7 X 10^° 


3 


0.3 


0.30 


2 


M33 


SA(s)cd 


56 


0.8 


1539 


3.3 X 10^° 


7 


0.36 


0.45 


3 


M51 


SAbc 


20 


9.7 


380 


1.1 X 10^^ 


17 


5.3 


3.1 


4 


M81 


SA(s)ab 


59 


3.3 


385 


1.2 X 10^^ 


8 


0.96 


0.89 


5 


M83 


SAB(s)c 


24 


3.7 


809 


3.4 X 10^^ 


18 


2.0 


1.6 


1 


NGC 253 


SAB(s)c 


78 


3.9 


2707 


1.3 X 10^^ 


35 


6.2 


6.3 




NGC 891 


SA(s)b? sp 


88 


9.6 


286 


8.1 X 10^^ 


6 


3.9 


3.3 




NGC 3628 


SAb pec sp 


89 


6.7 


247 


3.4 X 10^^ 


6 


2.0 


1.1 




NGC 4565 


SA(s)b? sp 


86 


12.5 


54 


2.6 X 10^^ 


7 


1.7 


1.3 




NGC 4631 


SB(s)d 


85 


7.5 


476 


8.2 X 10^^ 


6 


4.0 


2.1 




NGC 5775 


SAbc 


81 


26.7 


94 


2.1 X 10^2 


8 


10. 


7.3 




NGC 5907 


Sc 


87 


11.0 


72 


9.9 X 10^° 


4 


2.0 


1.3 




NGC 6946 


SAB(rs)cd 


38 


7.0 


457 


6.8 X 10^^ 


16 


3.3 


3.2 


4 



The columns are: source name, morphological type of a galaxy, i is the inclination (degrees), D is the distance of a galaxy (Mpc), 
S4.85 is the total flux density at 4.85 GHz, Bt is the total magnetic field strength, SFR is the star-formation rate estimated from the 
radio luminosity at 4.85 GHz, SFR (IR, Ha) is the star-formation rate estimated from IR or Ha emission, and the last column is the 
references to the SFR (Ha, IR) values. The references to SFR are: 1: Young et al. (1989), 2: Tabatabaei & Berkhuijsen (2010), 3: 
Verley et al. (2009), 4: Leroy et al. (2008), 5: Gordon et al. (2004). 



formation. For these galaxies, we integrated the observed 
radio intensities at 4.85 GH z and dete r mined their flux den- 
sities 54.85 as described in Istil et alJ (l2QQ9k (see Table [TJ. 



Dividing by the area of integration gives the mean total in- 
tensities /4.85 (mJy/beam) for these sources. We estimated 
the non-thermal intensities /n4.85 accounting for a mean 
thermal fraction of 20%. Further, from the 4.85 GHz flux 
density values 54.85 we determined the radio luminosities 
L4.85 and extrapolated these values to 1.4 GHz assuming a 
total spectral index of a = 0.7. The values for L1.4 are 
also summarized in Table [T] We used Eq. (l22l) to calculate 
the values for the SFR(radio) from the radio luminosity at 
1.4 GHz L1.4 (see Table [TJ. To check the reliability of esti- 
mated SFR(radio), we compared th em with the SFR va lues 
determined from the IR-emissio n dVoung et al.l 1989|) for 
these galaxies, adopting Eq. (3) o f lKennicuttl (ll998S or more 
recent values from the literature for special galaxies. The 
values determined from the radio emission and the other 
values correlate well for the 14 galaxies with SFR(Ha, IR) ex 
SFR(radio)^"^^^^'^^ with a correlation coefficient of 0.97 
(see Fig. [2]). Good agreement between independent mea- 
surements suggests that the SFR can be reliably estimated 
from the radio luminosity at 1.4 GHz (Eq. (|22]) ). 

Using the /n4.85 (see Table[T]) and the pathlength through 
the emitting medium (assumed to be the diameter of the 
galaxy for edge-on galaxies), we estimated the total mag- 
netic field strength 5tot under the assumption of the equipar- 
tition between the tot al energy of cosmi c rays and that of 
the magnetic field (Bec k & Krausel[2QQ5l. their Eq. (3), us- 
ing = 4.85 GHz and a synchrotron spectral index of an = 
0.8). The knowledge of the SFR of a galaxy and its effective 
diameter (or the pathlength through the emitting medium) 



allows the 5tot to be estimated, and vise versa, the SFR can 
be derived from the ^tot and effective radius of a galaxy. In 
turn, independent measurements of a SFR and 5tot open a 
possibility to calculate the effective diameter of unresolved 
sources in deep radio surveys. 

In the next section, we present simulations for one spe- 
cific SFR because a variation of the SFR simply scales with 
Btot in the case of an unchanged size of the galaxy during 
the evolution, as was assumed in this work. 

5 Modeling the regular magnetic fields and 
radio emission properties 

We perform modeling for a galaxy of the radius of 10 kpc in- 
clined at an angle of 60°, and having a constant SFR=10 Mq yr~^ 
which corresponds to the radio luminosity L1.4 = 1.8 x 
10^2 W Hz-i at 1.4 GHz (Eq. The strength of the to- 
tal magnetic field at the age t > teqi = 1.2 Gyr is estimated 
to be Btot (t) = 6.9 /iG, and the amplitudes of turbulent and 
regular fields are b = 6.2 /iG and B{t) = 3.1 /iG (Eqs. (O 
and (E])). 

We start the modeling of magnetic spots at the age of 
0.8 Gyr (Sect. l3.2l) . We assume that the formation of n = 15 
magnetic spots of a scale length of 1 kpc in the disk (with 
a scale radius 8.3 kpc and scale height 0.5 kpc) are in place 
at ^ 0.8 Gyr after the disk formation. The amplitudes of 
the regular field of magnetic spots are distributed randomly 
around the mean coherent magnetic field strength B{t = 
0.8) = 0.6 /iG, within the range from 0.7B{t = 0.8) to 
1.3B{t = 0.8), and coherent fields in the spots are oriented 
within the pitch angles 20° ± 5° or tt + (20° ± 5°) with 
equal probability. 
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Fig. 3 Simualtions in the framework of the SKA Design Studies (SKADS): the evolution of regular magnetic fields in 
the disk of a galaxy seen face-on (the frame size is 20kpcx20kpc). The amplitude and ordering scale of the regular fields 
are estimated at 0.8 Gyr after the disk formation epoch (~ 1.3 /iG and 1 kpc), 1.3 Gyr (~ 3.1 /iG and 1.6 kpc), 2 Gyr (3.1 
/iG and 2.5 kpc), 3 Gyr (- 1.5 /iG and 7 kpc), 5 Gyr (- 3.1 /iG and 6 kpc), and after 13 Gyr (- 3.1 /iG and 15.4 kpc). The 
color bar represents the amplitude of the regular field in units of /iG. 



Modeling the regular fields for an undisturbed (with- 
out a major merger) Milky- Way type galaxy (Fig. O show 
that the development of spotty structures with a coherence 
length of ~ 1.5 kpc is clearly visible during ^ 1.3 Gyr af- 
ter the disk formation. During the second Gyr the spotty 
patterns develop to prolonged structures of magnetic fields 
with a coherence length of 2.5 kpc with few reversals of 
the field. After ~ 5 Gyr the reversals disappear and the co- 
herence scale of the regular field reaches ~ 6 kpc. The fully 
coherent field is developed in ~ 5 Gyr while the asymmetry 
of the magnetic field strength remains longer. 

Simulations of the total intensity, polarization and Fara- 
day depth (see Appendix A for a mathematical description 
of equations) in the rest frame of a galaxy at 5 GHz af- 
ter 0.8 Gyr, 2 Gyr, 3 Gyr, and 13 Gyr of disk formation are 
shown in Fig. |4l The symmetric structure in total intensity 
at the age of 0.8 Gyr (Fig. IH first row) is due to strong 
turbulent field compared to the regular field (see Table |2). 
Smooth and weak polarized intensity and small Faraday ro- 



tation with no large-scale patterns are the result of a small 
ordering scale and small amplitude of the regular field at this 
age. At the age of 2 Gyr, the regular field has a coherence 
scale of ~ 2.5 kpc and few reversals, reaches the equipar- 
tition level and becomes comparable with the strength of 
the turbulent field (B ^ b/2). These changes are mani- 
fested by appearance of polarization patterns and asymmet- 
ric, inhomogeneous RM structures (Fig.lH second row). Af- 
ter 13 Gyr, the regular field has no reversals and is coherent 
at the size of a galaxy (Fig.|4l third row). The polarization 
map looks symmetric with notable depolarization along the 
major axis, and the RM structure is smooth and asymmetric. 
Note that the intensity units in Fig. |4l have to be scaled via 
SFR (or total radio luminosity; see Sect.©. 

Simulations in the rest frame of the star-forming galaxy 
at 150 MHz are shown in Fig. [5j The main difference be- 
tween the 5 GHz and 150 MHz maps is the effect of Fara- 
day depolarization which becomes stronger at low radio fre- 
quencies. Depolarization along the major axis of a galaxy 
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Table 2 The amplitudes of regular and turbulent fields, ra- 
dial and azimuthal coherence scales simulated for different 
ages of a disk galaxy with inclination of 60°, scale radius of 
R = lOkpc, and SFR^IOMq yr-^ 





B 


b 




let) 






(Gyr) 


(/iG) 


(/iG) 


(kpc) 


(kpc) 




(MHz) 





0.02 


6.2 


0.1 


0.1 


10 


1650 


0.8 


0.6 


6.2 


1.0 


2.9 


4.7 


855 


1.3 


3.1 


6.2 


1.6 


4.5 


3.6 


690 


2 


3.1 


6.2 


2.5 


6.8 


2.5 


525 


3 


3.1 


6.2 


3.6 


10. 


1.9 


435 


5 


3.1 


6.2 


6.0 


16.5 


1.1 


315 


13 


3.1 


6.2 


15.4 


42.4 


0.0 


150 



t is the age of a galaxy after the disk formation epoch assumed 
to occur at z = 10. 

z is the redshift of a galaxy calculated for a flat cosmology 
model (Qa + Qrn — 1) with Q.rn — 0.3 and Hq — 70 
kms~^ Mpc~^. 

z^z is the rest frame frequency which is detected at 150 MHz 
(A = 2 m) in the frame of the observer. 



is stronger at large ages when the regular field becomes 
stronger and more ordered. The depolarization shapes the 
elongated polarized structure aligned near the minor axis of 
a galaxy (Fig. [5]) which is clearly visible at 13 Gyr after the 
disk formation. 

We also simulate the total intensity, polarization and Fara- 
day depth in the rest frame of the observer at i^obs = 150 MHz 
(Fig. [6]) taking into account the frequency shift of a galaxy 
at a redshift z (see Appendix A for details). For more distant 
galaxies the synchrotron intensity, depolarization effect, and 
the regular magnetic field become weaker. The latter two 
effects lead to the small "observed" Faraday rotations at 
high redshifts (Fig. [6]). The observed intensity at 150 MHz 
is emitted at high frequencies at high redshifts (Table O. 

6 Discussion and conclusions 

The evolution of large-scale magnetic fields in galaxies has 
not been investigated in detail so far. Dynamo theory pre- 
dicts the generation of large-scale coherent field patterns 
("modes"), but the timescale of this process is comparable 
to that of the galaxy age (Sect. 2.2). Many galaxies are ex- 
pected not to host fully coherent fields at the present epoch, 
especially those which suffered from major mergers or in- 
ter actions with other galax ies. 

Gissinger et al.l ( 2009k presented the first global numer- 



namo action ( von Linden et al.l [199 8*) revealed large-scale 
field structures which are partly oriented along the spiral 
arms, but also features of regular fields in the interarm re- 
gions which rapidly change their structure within one rota- 
tion of the galaxy and hence do not lead to a stable pattern. 
Presently it is unknown whether feedback of the field onto 
the gas flow may stabilize the field pattern. 

The first global galactic- scale MHD simulation of a CR- 
driven dynamo wit h lOOpc resolution was performed by 
iHanasz et al in which the cosmic rays are produced 

in randomly occurring supernova remnants. The timescales 
for amplification and ordering of the magnetic filed in their 
m odel are consistent with the dynamo timescales estimtated 
bv lArshakian et al.l (l2009k . This gives some confidence that 
the dynamo model correctly describes the overall properties 
of field evolution, especially as the simulations by lHanasz et al 



ical simulation of a dynamo including superbubbles gen- 
erated by SNRs, but without explicit inclusion of the alpha- 
effect (no mean-field assumption). The generated spiral field 
is strong and of quadrupolar symmetry, confirming the mean- 
field results, and it resembles the "spotty" field injection. 

Kinematical simulations of the field evolution in spiral 
galaxies without feedback to the gas flow, but excluding dy- 



(l2009h also show "spotty" magnetic structures during the 
evolution from randomly spread seed fields to large-scale 
spiral structures. Even the ratio h/R regulates itself in the 
CR-MHD simulations to a value near 0.1 which we as- 
sumed in the dynamo models (Hanasz, private communi- 
cation). 

Spiral arms and self gravity of the gas need to be in- 
cluded in future improvements of the model. The inclusion 
of shock fronts from supernova remnants will need higher 
spatial resolution. Furthermore, the effects of the SFR on 
the strength of the amplification and saturation of the reg- 
ular magnetic field has still to be included into the model. 
As the SFR decreases with decreasing redshift and depends 
on merger rate, this may also influence the evolution of the 
magnetic field in evolving galaxies. 

Predictions of the model Comparison of our predictions 
with observations has to be restricted to spiral galaxies with 
thin disks and to high-frequency observations. Also, com- 
parison is still restricted to nearby galaxies which can be re- 
solved with the limited sensitivity of present-day telescopes. 
Among the existing observations, we could find spiral galax- 
ies with well-developed large-scale spiral ( axisym metric) 

magn etic field structures, like in M31 (lBeckll982: B erkhuiisen et al 
20031) . IC 342 (lKrauseetal.ll989h and NGC 6946 (lBeckl2007h . 
as visible in polarized intensity and Faraday rotation similar 
to the last row in Fig.Hl All three galaxies are large with an 
ex ponential scale radius of about 5 kpc. Hence, according 



ex ponential scale radius oi about ^ Kpc. Hence, accoraing 
to lArshakian et al. I (l2009k it should have taken about 9 Gyr 



for these galaxies to coherently order the regular field struc- 
ture over the entire disk (twice scale radius). This indicates 
that they did not suffer any major merger during this period. 
Indeed, NGC 6946 and IC 342 are isolated galaxies without 
any known companion and without any signs of major tidal 
distortions. 

Observations of M 5 1 showed, however, that even if a 
large-scale spiral magnetic field is present in the galaxy's 
disk and apparent in pola rized intensity, it i s hidden in the 
map of Faraday rotation (F letcher et al and looks much 

more patchy than that in row 3 of Fig. |4] This may be due to 
an additional magnetic field component, probably an isotropic 
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Fig. 4 The rest frame simulations of the total intensity {left panel), polarization (middle panel), and Faraday rotation 
(right panel) at 5 GHz for a galaxy with an inclination angle of 60°, turbulent (6.2 jiG) and regular (from 1.3 jiG to 3.1 /iG, 
see Table [21) magnetic fields, and star-formation rate of 10 Mq yr~^ are shown for 0.8 Gyr, 2 Gyr, 3 Gyr, and 13 Gyr after 
disk formation. The frame units are given in kpc. The color bars in the first and second columns (total and polarized 
intensity) are given given in arbitrary units. The color bar of the third column (Faraday rotation measure) is given in units 
of radm~^. 



turbulent magnetic field with scales of the longer axis be- 
tween 400 pc to 1 kp c. Such a field has been found in M 51 
([Fletcher et al."2010') and in the barred galaxy NGC 1097 
dBeck et al. 2005), but was not included into our simula- 
tions shown in Fig. |4] and Fig. [5j This component may be 
produced by strong shearing gas motions or compressions. 
It may also form loops perpendicular to the galaxy's disk 
due to the Parker instability. 

A prediction of our model is that interacting and merg- 
ing galaxies should reveal complicated field patterns. The 
nearby and best studies case of strongly interacting/merging 
galaxies is the Antennae galaxy pair N GC 4038/39 observed 
in total and polarized radio emission by' Chyzy & Beckl (|2004|) 
A huge polarized ridge in the northeast is probably the result 
of shearing flows. Another highly polarized region between 
the two galaxies marks the location of strongly compressed 
fields. The spiral field pattern of the galaxy in the northwest 
is still clearly visible. This galaxy is still in the phase be- 
fore total field disruption and cannot be compared with our 
simulations. 



A survey of polarized emission fro m g alaxies in the Vir- 
^o clu ster by lWezgowiec et all ( 2007 ) and lWezgowiec et al 



revealed signs of compression or shear in almost all 
galaxies, but the Faraday rotation data are still of insuffi- 
cient quality to detect any delay in the formation of coherent 
fields. 

Another prediction from our model, that small galaxies 
build up their large-scale field faster than large galaxies, also 
cannot be observationally confirmed with the existing data 
from nearby spiral galaxies. We need a sample of galaxies 
without massive density waves and without bars where the 
dynamo-generated field dominates over anisotropic fields. 
Furthermore, galaxies with strong signs of tidal interactions 
have to be excluded. A comparison with observations of 
dwarf and irregular galaxies fails because they do not fulfill 
the "thin-disk criteria" (h/R < 0.1) discussed in Sect. 13.51 

Forthcoming observations of a large sample of spiral 
galaxies with high resolution and high sensitivity and of dis- 
tant galaxies in early stages of the evolution with the SKA 
and its precursors should be able to test our scenario (.Beckl 
[2009h . 
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Limitations of the model. Lacking a realistic MHD model 
of evolving galaxies, we tried in this paper to describe the 
evolution of large-scale fields by simple interpolation be- 
tween the stages about which we believe to have some knowl- 
edge: the epoch of field seeding and the present epoch. This 
approach provided templates for radio maps in total and po- 
larized intensity and in Faraday rotation measures, as part 
of the simulations for the SKA Design Studies (SKADS). 

In Arshakian et al. (2009), we assumed that the "quasi- 
spherical" mean-field dynamo amplified the regular field 
and increased the coherence scale in radial direction, while 
the "disk" mean-field dynamo was effective in thin-disk galax- 
ies. We note that the simulations of the "disk" mean-field 
dynamo described in the previous section cannot be per- 
formed for dwarf and irregular galaxies. 

Another limitation of simulations comes from the fact 
that the strength of the cosmic microwave background en- 
ergy density increases with increasing redshift, so that the 
synchrotron emission of star-forming galaxies will be quick- 
ly suppressed by inverse Compton losses off the cosmic mi- 
crowave background at redshifts z > 3 (Murphyl l2009k . In 
this paper, we do not account for this effect, but our sim- 
ulations should be realistic for star-forming galaxies up to 
z^3. 



High SFR causes high velocity turbulence of the ion- 
ized gas in starburst galaxies, e.g. by a major merger of gas- 
rich galaxies. If SFR is higher than or equal to 20 yr~^, 
the mean-field dynamo will be suppressed (Arsh akian et al. 
12009) . Hence, our simulations are valid only for disk galax- 
ies with SFR ^ 20 Mq yr'^. 

In reality, the magnetic field evolution in first galax- 
ies can be affected by various effects which are not ac- 
cou nted in our simple model. In particular, g alactic winds 
(e^Moss et al.l2010h . turbulent pumping re.g. lBrandenbur^ et al. 

and the multiscale nature of the interstellar medium 
(iMantere et al can substantially affect the dynamo 

and result in a phenomenology which observably differs 
from predictions of our scenario. Future observations with 
SKA and its forthcoming pathfinders will be crucial to de- 
tect the difference with our simplified model and prove the 
importance of additional effects. 

Furthermore, our simulations assume that the field struc- 
ture in the halo is the same as in the disk and do not take 
into account the X-shaped halo fields observed around edge- 
on galaxies ( Krause 2009). When observing mildly inclined 
galaxies at low frequencies, the polarized emission from the 
disk is mostly depolarized and the halo may dominate. In 
this case, our results for polarized and Faraday rotation im- 
ages shown in Fig.[5]may change significantly. This can be 
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tested by sensitive observations of edge-on galaxies at suf- 
ficiently high frequencies to avoid Faraday depolarization. 

In summary, we link the SFR and amplitudes of the reg- 
ular and turbulent fields in disk galaxies, further develop 
an evolutionary model for magnetic fields in isolated star- 
forming disk galaxies and performed modeling of the evolv- 
ing galaxy. The amplitude of the initial "spotty" structure of 
regular fields (generated at the epoch of disk formation) is 
amplified by means of mean-field dynamo to the equiparti- 
tion level in a few Gyr and remains at this level up to the 
present time. The ordering and azimuthal scales of initial 
magnetic spots increase with the age of a galaxy stretching 
the initial field in both directions. For the proposed model, 
we simulate the total intensity, polarization and Faraday ro- 
tation for the MW-type galaxy at different frequencies (from 
150 MHz to 5 GHz) observable with the SKA and its precur- 
sors. A number of predictions such as the patchy structure 
and field reversals of regular fields, patchy Faraday rota- 
tions in galaxies younger than few Gyrs, smaller ordering 
scale in young galaxies, polarization patterns at 5 GHz and 
asymmetric structure at 150 MHz in older galaxies because 
of stronger depolarization, and the complicated field pat- 
terns in interacting and merging (minor and major) galaxies 
can be tested with the future radio telescopes. 
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A Equations for a total and polarized 
intensities, and Faraday rotation 

The polarized signal passes a distance w' in the direction of the 
line-of- sight through a magneto-ionic medium of the galaxy in- 
clined at an angle i. Then the Faraday rotation (in units of rad 
m~^), polarization angle, and total intensity in the coordinate sys- 
tem inclined at an angle i {x,y' = y cos i,w' = w cos i) are 

B\\{r)ne{r)ds, (Al) 

-oo 

x{x, y', w', A, z) = xo{x, y' , w')+ 

and 

I{x,y',\,z) = {y^Y J ne(r)|Bx(r)|^+°ds, (A3) 

where B is the magnetic field strength measured in /xG, rie in 
cm~^, ds and w' in kpc, xo is the intrinsic magnetic field orienta- 
tion (from the magnetic field model), A is the observed wavelength 
(in meters), z is the redshift of a galaxy, a is the synchrotron spec- 
tral index (a = 0.8). The equations for Stokes parameters, Q and 
U, are given by 

Q{x,y\\z) = 

p\yV~z) J Mr)\B^{r)\^^'' cos{2x)ds, (A4) 
U{x,y',\z) = 

P^[yT~z) / Mr)\B^{r)\'^'^s\n{2x)ds, (A5) 

— oo 

where pi — 0.73 is the intrinsic maximum polarization for syn- 
chrotron emission (for a — 0.8) in regular magnetic fields. Then, 
the "observed" polarization intensity 

P/(x, 2/^ A) = ^Q{x,y',\Y + U{x,y',\)\ (A6) 

Xo.s{x^y\X)-\^rci.n (A7) 
and the "observed" RM is 

iiMobs(x, y', X) = (A8) 
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